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Section  1 
INTRODUCTION 


1.1  THE  BASIC  PROBLEM 

J  The  problem  addressed  by  this  exercise  is  the 
extraction  of  a  sequence  of  points  C^ranges*^  along  specified 
tracks  (generally  radials  emanating  from  a  given  point)  which 
intersect  digitized  contour  lines  (of  e.g.,  bathymetry), 
along  with  the  associated  contour  values.  The  points  are 
then  sorted  to  provide  the  contour  values  as  a  function  of 
range  along  the  sequence  of  radials. 

The  digitized  contours  and  radials  are  specified  as 
a  sequence  of  (x,y)  [longitude,  latitude]  pairs.  A  Cartesian 
geometry  is  assumed.  Since  distances  are  computed  in  the 
standard  Cartesian  way,  it  is  totally  irrelevant  whether  the 
radials  and  contours  are  specified  in  terms  of  an  absolute 
origin  or  in  terms  of  a  local  relative  origin.  This  latter 
feature  of  course  has  advantages  when  the  absolute  coordi¬ 
nates  are  sensitive.  In  addition,  the  method  makes  no  as¬ 
sumption  about  whether  the  contours  are  closed  or  open  in  the 
"universe"  under  consideration. 


Section  2 

MATHEMATICAL  METHOD 


The  radials  under  consideration  in  this  problem  are 
considered  to  be  simple  one-component  line  segments.  The 
digitized  contours  are  considered  to  be  a  sequence  of  line 
segments,  which  may  or  may  not  be  closed.  The  mathematical 
solution  of  the  intersection  of  a  line  (radial)  with  a  se¬ 
quence  of  segments  (contour)  depends  on  an  interesting  pa¬ 
rameterization  of  a  line  segment.  Given  a  line  segment 
defined  by  its  endpoints  (X'j  and  y^)  and  {%2  'i2^  • 

we  choose  to  express  the  equation  of  the  line  in  the 
following  parametric  form. 


t^  >  1,  the  point  (x,y)  is  somewhere  on  the  extension  of  the 
original  segment  but  external  to  it. 

We  will  now  parameterize  a  second  line  segment, 
defined  by  its  endpoints  (X3»y3)  and  (x^jy^)  in  a 
similar  fashion,  using  the  parameter  t2* 


^  ^2  *  ^*4  "  *3^ 


y  ■  y^  +  t2  *  (y^  -  y3)  •••  the  second  line  segment 
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Now  equate  the  two  values  of  x  and  the  two  values  of  y  and 
solve  for  the  values  of  ti  and  t2*  We  have  the  following 

-(x^-x^ ) (y^-y3)  +  (y^-y^ ) ( x^-x^ ) 

"  -(x,-x, ) (y .-y,)  +  {y,-yi ) (x.-x,) 


(X2-X3  )  {y3-y^ )  -  (y2-y2)  ) 

^2  "  -(X2-X^ )  (y^-y3)  +  {y2-y2) 


Note  that  the  denominators  are  equal.  If  the  common  denom¬ 
inator  is  zero,  then  the  two  line  segments  are  parallel  and 
do  not  intersect.  See  case  a)  below.  If  both  t-i  and  t2 
are  greater  than  zero  and  less  than  1,  then  the  line  segments 
intersect.  See  case  b)  below.  If  either  t2  or  t2  are 
less  than  zero  or  greater  than  1,  then  the  two  segments  do 
not  intersect  along  the  extent  of  their  definition.  See  case 
c)  below  for  some  examples. 


(a)  ^  (b)  -  (c) 

We  simply  apply  this  test,  checking  a  given  radial  in  turn 
against  each  line  segment  making  up  each  contour.  Remember 
that  simple  Cartesian  geometry  is  used.  Also  note  that  there 
is  no  requirement  that  the  contours  be  closed. 


DATA  BASE  FORMAT 


The  file  of  digitized  contours  will  necessarily  be 
supplied  from  an  existing  data  base.  In  the  present  computer 
implementation,  the  following  format  is  assumed.  There  are 
four  records  in  the  formatted  file  for  each  contour.  Their 
content  is  as  follows. 

Record  1.  FORMAT( IX, 13)  IVAL  ...  The  contour  value 

Record  2.  FORMAT (AlO)  X  ...  A  dummy  value 

Record  3.  FORMAT( 4X,  14)  NPTS  ...  The  k.umber  of 

points  in  the 
contour 

Record  4.  FORMAT ( 2P10. 2)  (ALON(I),  The  LONGITUDE-LAT- 

ALAT( I) ,1=1 ,NPTS)  ITUDE  pairs  making 
...  up  the  contour 

When  changing  or  modifying  this  code  for  another  data  base, 
it  will  be  necessary  to  modify  the  four  statements  in  the 
present  code.  The  changes  may  include  all  or  some  of  the 
following. 

a.  Modification  of  format  to  accept  binary  input. 

b.  A  change  in  the  order  of  the  parameters. 

c.  The  number  of  READ  statements  necessary  to  ac¬ 
complish  the  task. 

The  input  section  is  thus  seen  to  contain  no  serious  restric¬ 
tions  on  the  format  of  the  data  base.  Note  though  that  the 
contours  are  presented  totally  one  after  the  other,  each  is 
complete  in  itself.  No  other  order  of  contour  sequencing  is 
expected  or  required  by  the  code. 


COMPUTER  IMPLEMENTATION 


The  code  as  presented  in  the  appendix  is  straight¬ 
forward.  The  following  will  serve  as  a  simple  narrative  of 
the  computation. 


a.  (Fortran  label  10)  Get  the  contour  from  a  disc 
file.  As  mentioned  above,  this  segment  is 
likely  to  change  in  another  case.  Skip  to  the 
end  of  the  computation  section  (Fortran  label 
80)  if  there  are  no  more  contours  to  be 
considered. 


b.  Determine  contour  limits  in  spatial  terms. 


c.  DO-loop  on  radials  ...  DO  50  NR  =  1,  NRADS. 


d.  Determine  physical  limits  of  the  given  radial; 
skip  consideration  of  this  radial  if  there  is 
no  possible  intersection. 


e.  DO-loop  on  contour  segments  ...  DO  40  1  =  1, 
NMI,  where  NMI  is  defined  to  be  one  less  than 
the  number  of  points  in  the  contour  under 
consideration. 


f.  Check  for  possible  intersection  of  given  radial 
with  current  segment  of  contour.  Store 
relevant  results  if  so. 


g.  End  DO-loop  on  points  in  this  contour  (label 
40)  and  on  the  number  of  radials  under 
consideration  in  this  problem  (here,  fixed  at 
five) . 


h.  (Fortran  label  80)  Sort  intersections  for  each 
radial  into  increasing  order  of  range  and  then 
return. 


VARIABLE  Dictionary 


ALON(IOOO) 


ALAT(IOOO) 


NPTS 


X1S{5) , 
Y1S( 5) 


X2S(5) , 
Y2S(5) 


longitude  values  for  a  given  contour 


latitude  values  for  a  given  contour 


number  of  points  in  current  contour 


starting  longitude  and  latitude 
values  for  each  of  five  radials 


final  point  {longitude  and  latitude) 
for  each  of  five  contours 


SAVRAD( 300 , 5 )  for  each  of  300  total  radial-contour 
intersection  in  the  problem  the  fol¬ 
lowing  : 


longitude  value  (-,1) 
latitude  value  (-,2) 
range  (nm)  (-f3) 

radial  ordinal  (-,4)  (in  the  range 


1-5) 


contour  value  {-,5) 


4-2 


Section  5 

POSSIBLE  EXTENSIONS  TO  THE  CODE 

There  are  always  many  possible  extensions  or 
modifications  to  a  computer  code  as  new  requirements  arise. 
The  following  is  a  list  of  some  of  the  more  likely  ones  for 

this  code. 

a. 

The  number  and  dimension  of  the  radials 
permitted  should  be  input  to  the  code  in  its 
role  as  either  a  main  program  or  as  a  sub¬ 
routine. 

b. 

The  limit  of  300  total  radial-contour 

intersections  is  likely  to  prove  restrictive. 

c. 

(In  conjunction  with  b)  immediately  above)  A 
cutoff  on  the  maximum  range  of  a  radial-contour 
intersection  may  be  desirable  in  some  cases. 

d. 

The  code  segment  towards  the  beginning 
governing  the  specific  form  of  the  digitized 
contour  database  will  almost  certainly  be 

changed . 
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